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Abstract 

Black hole solntions of type IIB snpergravity have been previonsly con- 
strncted that describe the N=4 snpersymmetric Yang-Mills plasma with 
an anisotropic spatial deformation. The zero temperature limit of these 
black holes approach a Lifshitz-like scaling solution in the infrared. We 
show that these black holes become unstable at low temperature and 
we construct a new class of black hole solutions which are thermody¬ 
namically preferred. The phase transition is third order and incorpo¬ 
rates a spontaneous breaking of the SO{6) global symmetry down to 
S'0(4) X SO{2). The critical exponents for the phase transition are given 
by = (—1,1,1,2) which differ from the standard mean-held 

exponents usually seen in holography. At low temperatures the black 
holes approach a novel kind of scaling behaviour in the far IR with spa¬ 
tial anisotropy and hyperscaling violation. We show that the new ground 
states are thermal insulators in the direction of the anisotropy. 



1 Introduction 


Strongly coupled conformal field theories that have been deformed by spatially de¬ 
pendent sources in fiat spacetime can be studied by constructing novel black hole 
solutions using the AdS/CFT correspondence (e.g. Such studies are interest¬ 

ing for a number of reasons. Without such sources the translational symmetry im¬ 
plies that momentum cannot dissipate and this leads to non-physical delta-function 
responses in, for example, the thermal and electric conductivity of the system. Spa¬ 
tially dependent sources provides a mechanism for momentum to dissipate and this 
leads to hnite DC responses. The spatially dependent sources also provide a useful 
tool to search for novel holographic ground states which can appear in the far IR; 
insulators, coherent metals and incoherent metals have been realised in this way, as 
well as transitions between them . A more specihc motivation derives 

from the properties of the quark gluon plasma observed in heavy ion collisions. In 
particular, the plasma appears to have regimes where it is described by a strongly 


coupled and spatially anisotropic fluid 26,27 


An interesting framework for analysing spatial anisotropy in N=4 super Yang- 
Mills (SYM) theory was initiated in and then further developed in [^[^. Specif¬ 
ically, black hole solutions of type IIB supergravity were constructed in that 
asymptotically approach AdS^ at the UV boundary with the type IIB axion having 
a linear dependence on one of the three spatial coordinates. The linear axion source 
is associated with a distribution of D7-branes that intersect D3-branes in two of the 
spatial directions and is smeared in the third. At low temperatures these black holes 
approach a T = 0 solution, constructed in [^, which becomes a Lifshitz-like scaling 
solution in the far IR. 

It is natural to interpret this scaling solution as the T = 0 ground state of the 
anisotropically deformed N=4 SYM theory. Here, however, we will show that the 
black holes of are unstable at low temperatures and there is a phase transition 
which spontaneously breaks the global 50(6) symmetry down to 50(4) x 50(2). The 
origin of this instability was already noticed in (^. In particular, by analysing the 
Kaluza-Klein spectrum of the five-sphere it was found that there are scalar modes, 
transforming in the 20' of 50(6), which saturate the BF bound in AdS^ background 
but violate an analogous bound in the Lifshitz-like background. This suggests that 
the Lifshitz-like scaling solution is unstable. It is natural to suspect that such an 
instability is also present for the T = 0 solution that interpolates between AdS^ in 
the UV and the Lifshitz-like solution in the IR. By continuity one then expects that 
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the finite temperature black hole solutions should become unstable at some finite 
temperature. 

Here we will show that these expectations are realised. We find that the black 
holes become unstable at some critical temperature with two new branches of 
black holes appearing. One of these branches seems to be a non-physical branch of 
“exotic hairy black holes” 28,29 . Specifically, it seems that these black holes only 


exist for temperatures T >Tc and are never thermodynamically preferred. The other 
branch exists for T < and is thermodynamically preferred. The phase transition 
is continuous and, somewhat surprisingly, third order. Furthermore, we calculate 
the critical exponents of the phase transition finding {a, f3, 7 , 5) = (— 1 , 1 , 1 , 2 ) rather 
than the standard mean-field values {a,(3,'y,6) = (0,1/2,1,3) associated with most 
holographic phase transitions. The critical exponents we find have been previously 
realised in bottom-up models of holographic superconductors with a Lagrangian con¬ 


taining cubic terms in the modulus of the complex scalar field 30 33 . While such 
terms are a bit unnatural for a complex scalar field they are natural for a neutral 
scalar field, provided the potential does not have any discrete symmetry. The critical 
exponents that we find can, in a certain sense, be realised by a Landau-Ginzburg 
(LG) model with a scalar order parameter and cubic term in the free energy. While 
such cubic terms in LG models are associated with first order phase transitions, our 
holographic transition appears to be continuous and third order. 

We construct our new solutions using a consistent KK truncation of type IIB 
supergravity on that keeps the D = 5 metric coupled to the axion and dilaton, as 
in , and in addition keeps an extra single neutral scalar field. Any solution of the 
D = 5 theory gives rise to an exact solution of type IIB supergravity. The potential 
for the neutral scalar field in the D = 5 theory does not have any discrete symmetry 
and this is associated with the non-standard values of the critical exponents that we 
just discussed. The fact that the phase transition spontaneously breaks SO(6) to 
S'0(4) X SO{2) is only apparent after uplifting to type IIB. 

We construct the new branch of black hole solutions down to low temperatures and 
elucidate the T = 0 behaviour. Similar to and unlike many holographic studies, 
the IR part of the geometry at T = 0 does not approach a scaling solution of the 
equations of motion, but instead approaches the leading terms of an expansion, which 
eventually approaches AdS^ in the UV. The leading terms of this IR expansion are 


similar to the hyperscaling violation solutions 34 36 but with anisotropic scaling 
in the spatial direction rather than the time direction. We calculate the thermal 
conductivity of the black holes, essentially importing the results of [16]. We show 
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that the scaling behaviour implies that at low temperatures the system is a thermal 
insulator with k ~ 2"io/3_ 

The remainder of the paper is organised as follows. In section we present 
the D = 5 top-down model and show how it is arises from a KK reduction of type 
IIB supergravity on S^. The construction of the black holes and a study of their 
properties is contained in section We conclude in section and we have two 
appendices. Appendix A contains some technical results concerning a Smarr formula, 
while appendix B discussed how the critical exponents that we hnd can be extracted, 
in a certain sense, from a Landau-Ginzburg type analysis. 


2 The top-down model 


We will consider a D = 5 gravity theory coupled to three scalar helds, the axion and 
dilaton, 0 and y, as in [^[^, and an additional scalar X. The bulk action is given by: 

S=j {r - ?,X-\dXf + 4(X2 + 2X-1) - , (2.1) 

where for simplicity of presentation we have set IbvrG = 1. The corresponding 
equations of motion are given by: 


= 0 , 

R,, = 3X-^d,Xd,X - + \e^^d,xduX • ( 2 - 2 ) 


This top-down model arises as a consistent truncation of the Kaluza-Klein (KK) 
reduction of type IIB supergravity on a hve-sphere. That is, any solution to the 


equations of motion (2.2) gives rise to an exact solution of type IIB supergravity with 
D = 10 metric and self-dual hve-form given by: 


dsio = + XA^Ad ^2 + X^A-^/^ ^ ^-l/ 2 j^-l ^Qg 2 ^ 

F( 5 ) = 2 t/vol 5 -I- 3 sin^ cos^X”^ *5 dX A d^ 

+ A“^sin^ cos^ ^{2Ud^ — 3 sin ^ cos^X~‘^dX) A dr A V0I3 , (2.3) 


where dsf and V0I5 are the D = 5 metric and volume form, respectively, dfl^ and V0I3 
are the metric and volume form on a round three-sphere, respectively, and 

A = X-^sin^^ + Xcos^C, K = X^ cos^^ + X”^ sin^^ + X'G (2.4) 
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The D = 10 dilaton and axion are the same as the D = 5 scalar helds 0 and X) 
respectively, and the H = 10 three-forms are both zero. When X ^ 0 this class 
of H = 10 metric and hve-form has the S'0(6) symmetry of the round hve-sphere 
reduced to S'0(4) x SO (2), with the first factor acting on the round and the second 
acting on the circle parametrised by r. 


That this is a consistent truncation can be established using the results of 37 


Indeed it was shown in that there is a consistent KK truncation of type IIB 
supergravity on a hve-sphere to Romans D = 5 SU{2) x f/(l) gauged supergravity, 
whose bosonic helds consist of a H = 5 metric, a scalar X, SU{2) x 17(1) gauge-helds 
and two two-forms. This truncation can simply be extended to include the Zl = 10 
axion and dilaton and we can then truncate away the gauge-helds and the two-form 
to obtain our model. 

Notice that the unit radius AdS^ vacuum solution to the equations of motion 


(2.2) has X = 0, X = 1 with constant 0, and uplifts to the standard AdS^ x 


solution of type IIB. Around this vacuum solution, perturbations of the helds 0, x 
are massless and are associated with marginal operators in X = 4 SYM theory with 
scaling dimension A = 4. Perturbations of X have = —4, which saturates the 
BF bound, and is associated with an operator with dimension A = 2. This 
operator is part of a multiplet, transforming in the 20 ' of 50(6) which is dual to 
operators in N=4 SYM constructed from the six adjoint scalar helds 0^ of the form 
Tr(0-^0‘^) — trace. When acquires an expectation value spontaneously, as it will 
in our solutions, it breaks the 50(6) global R-symmetry down to 50(4) x 50(2). 

Notice that setting X = 1, which is a further consistent truncation, we have 


A = 1, 17 = 2, from (2.4), and we recover the 0 = 5 model that was studied in 


In particular the metric on the hve-sphere in (|2.3|) becomes the round metric. Thus, 

to include one extra scalar 


our model extends the top-down model studied in 
held, X, which saturates the BF bound. It is sometimes convenient to consider a 
canonically normalised scalar held, 0, instead of X, dehned by 


A' 


(2.5) 


in terms of which the action reads 

S= J d^x^/^ (^R - ^(00)^ + 4(e-27/v^ + 2e^/^) - ^(00)^ - ■ 

( 2 . 6 ) 
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2.1 Brief review of previous work 


The anisotropic solutions constructed in Si , with the axion linear in one of the 
spatial coordinates, all lie within the ansatz 


ds^ = 


e 2 '- 




dv? 

J'Bdt^ H—— + dx'^ + dy'^ + Tddz^ 


7 


X = az, 


<t> = <l>(u). 


(2.7) 


with trivial X-field, X = 1 (i.e. -0 = 0 in (2.5)). The functions X, B are functions of 
the radial coordinate u, and the function "H is taken to be "H = e~'^ which, remarkably, 
can be imposed consistent with the equations of motion. 

The black hole solutions constructed numerically in approach in the UV, 
located at u —)■ 0, a unit radius AdS^ with a linear axion deformation. As T —)■ 0 the 
black hole solutions approach a T = 0 domain wall solution whose IR limit approaches 
a fixed point solution with 


_ 49 fl2\^ 2 


B = X- 


= 


2 

11\5 _4 

12 ' “ " 


( 2 , 8 ) 


which was first found in After switching to a new radial coordinate u = r 
this solution can be written in the form 


ds^ = L^[^+r^ 


-dt^ + dx"^ + dy‘^) + 


^ = az^ 


(2.9) 


with X = 1, where the bars denote quantities that have been rescaled, and lA = 
11/12. This metric is manifestly invariant under the anisotropic Lifshitz-like scaling 
{t, X, y, z, r) —)■ (At, Ax, Xy, A^/^A, A“^r). 

Following we can study the properties of a massive scalar field, satisfying 


V^cr = m^cr, in the background (2.9). By considering solutions of the form and 
demanding that A± are real, we deduce that m? > —11/3. Since we would like to 
identify A± as scaling dimensions in a putative field theory dual to these Lifshitz- 


like solutions, this suggests that the anisotropic solution (2.9) will be unstable under 
perturbations by any massive field with < —11/3. If the Lifshitz solution is 
unstable we expect that the T = 0 domain wall solution itself will be unstable and 
hence, by continuity, that the finite temperature black hole solutions will be unstable 
up to some critical temperature T^. 

Thus, since the scalar field X in our model has = —4, we anticipate that 
the black hole solutions of will still describe the high temperature phase of 
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the system but will become unstable at leading to a phase transition. The critical 
temperature can be found by establishing the existence of a suitable zero-mode in the 
linearised fluctuations of the X-£eld about the numerically constructed black holes 
of [2j[^. We carried out this analysis but we will omit the details. Instead we will focus 
on the construction of the new branch of fully back reacted black holes that emerge 
at T = Tc and examine some of the physical properties of the new low-temperature 
phase. 


3 Construction of new anisotropic black holes 

3.1 Ansatz and equations of motion 

We extend the ansatz of by allowing for a non-trivial X-£eld and consider 

^ (—tFBdt^ + dx^ -|- dy^ + Tidz^ H—, 

\ tF J 

X = az, 0 = 0(m), X = X{u), (3.1) 


where X, B and 1-i are functions of u. The function Fi is associated with the anisotropy 
in the 2 ; direction that is sourced by the axion held. By combining the equation of 
motion for the dilaton with the Einstein equations, as in [T]-[^, we find that it is 
possible to choose the function FL to be related to the dilaton via 


n = e-^, 


(3.2) 


and we will do so in the sequeQ 

We now discuss the resulting equations of motion for this ansatz, following the 
approach of [^. The equation for the axion y in (2.2) is trivially satished. The X 
equation of motion implies that 


Uu'^XXX" + 3^ - bu^XX4>' + 2u^XX^ + Au^XX' - UuXX^X' 

- 12u^X {X'f + lQe-‘^^^X{X^ - 1) = 0 , (3.3) 


while the 0 equation of motion gives 

B' 

AuX(^" + 2uX—(t)' + ^uX'ct)' - buX (0')^ - = 0 . (3.4) 

^Our preliminary investigations into relaxing this condition did not reveal any other solutions of 
physical interest. 
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There are also four independent components of the Einstein equations arising from 


(2.2). By taking a suitable combination of one of these equations with the 0 equation 


of motion, in order to eliminate B' terms, we can arrive at an equation which can be 
algebraically solved for 




e 2^ 


4(0' + M0" 


6 2*^0 ^ (4m + + 160') + 


4e-i'^(l -X)2(2 + X)0' 
3X(0' + m0") 


(3.5) 


Next, by taking a suitable combination of two of the remaining three Einstein equa¬ 
tions, in order to eliminate J^'' terms, we arrive at an equation that we can solve for 
B'lB: 


B' _ (240' - 9m0'2 + 2Om0" 
~B ~ 24 + 1Om0' 


24mX'2 


X2 (12 + 5m0') ’ 


(3.6) 


and we observe that only the combination B'jB appears in (3.3) and (3.4). It is 


now possible to show that (3.3), (3.4), (3.5) and (3.6) imply that all of the Einstein 
equations are solved. To see this we can use (3.4) to solve for X' in terms of 0 and 


X and their derivatives, after using (3.5) and (3.6). Furthermore, comparing this 


equation with the expression for X' that can be obtained by differentiating (3.5), we 


obtain a third order equation for 0, which can be used instead of (3.4). One can 


then check that the remaining Einstein equations are satisfied. Observe that if we 
set X = 1 we recover the equations of motion given in [^. 

In summary, the equations of motion are equivalent to ( |3.3 )-(3.6) and are, effec¬ 
tively, second order in X, third order in 0 and hrst order in B^ with X algebraically 
specihed by 0, X and their hrst and second derivatives. Thus, a solution is specihed 
by six integration constants. 

We note that the ansatz and hence the equations of motion are invariant under 
the following two scaling symmetries 


M — )■ Am, {t, x,y, z) ^ x,y, z), a —)• A ^a; 

( ^ Ai, B -» A-’/^B; 


(3.7) 


where A is a constant. 
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3.2 The UV and IR expansions 


We now discuss the boundary conditions that we will impose on (3.3)-(3.6). In the 
UV, as M —)■ 0, we demand that the asymptotic behaviour is given by 

4 (121a^ + 1152134 + 2304(X2)2^ 




+ u 


4032 


- u \ogu— + 

D 


lla^M^ 4 4 

= 1 H-^4- Vu J^ 4 ^ + u log u—— + ... , 


24 

^ lla^M^ 4, 

B =1 -- h M^i34 — logM—— + 


24 

X = 1 + U^X2 — u 


12 

7a^ 

12 


^ 50^X2 
24 


+ 


(3.8) 


The solutions are asymptotically approaching AdS^ with an anisotropic deformation 
of the axion field in the z-direction with strength a. This UV expansion is specified 
by four parameters, X4,134, X2, whose physical interpretation will be discussed below, 
and a. It is important to observe that we have set a possible n^logti term in the 
expansion of X to zero, as this would correspond to sourcing the operator dual to 
X which we don’t want[^ We next note that the second scaling symmetry in (3.7) 
has been used to set the leading term in B to unity. We also observe that associated 
with the first scaling symmetry in (3.7) the UV expansion is preserved under the 
transformations u ^ Xu and 


.-4 


134 ^ A 
X4 ^ A-V4 
X2 —)■ A ^Xo 
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i34+-aW-MogA, 

^a^A"^ log A, 


(3.9) 


The presence of the log terms is associated with the fact that the linear axion defor¬ 
mation gives rise to a non-vanishing conformal anomaly, as discussed in [^[^. 

In the IR, we will assume that we have a regular black hole horizon located at 


2 


In section 3.5 


when we calculate critical exponents of the phase transition, we will briefly 


consider black holes with such a source for X. 
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u = Uh- We therefore will demand that as u ^ Uh we have 
12uhXha?e^ 


— (j^h — 




32 + 3a‘^e2ulXh + 
X = Xh{u - Uh) + ... , 


{u-Uh) + . 


„ „ , 2B,(i5aV^"-utXi - 266(.Y3 - 1)^ - 9ea^i^ulX,(2 + X^)) , 

o — Oh H-:--r "9 -[u — Uh) + ■■ ■ , 


X = Xh 


Uh (32 + 3Xe-^ulXh + leX^y 
lQXh{Xl - 1 ) 


{u-Uh) + ... . 


(3.10) 


Uh (32 + 3 a?e-^ulXh + lOX^) 

This IR expansion is specihed by four parameters, (l)h,Bh,Xh and Uh, with Xh hxed 


via 


e-^ (32 + iii^e-^ulXn + 16 xh 

TH = r(uH) = -i-r——-(3,11) 

llUhJ^h 

We have noted that the equations of motion are specihed by six integration con¬ 
stants. We have eight parameters appearing in the asymptotic expansion minus one 


for the remaining scaling symmetry in (3.7). We thus expect to hnd a one-parameter 


family of solutions which can be parametrised by the quantity T/a. We note that the 
presence of the conformal anomaly introduces an additional dynamical scale which 
we hold hxed to be unity throughout our discussion. 

3.3 Stress tensor and thermodynamics 

To calculate the free-energy and the stress tensor, we need to supplement the bulk 
action with boundary counter terms (e.g. 1^). We write 


Stotal Sfjulk T Set , 


(3.12) 


where Sbuik is the bulk action given in (2.1) (or (2.6)) and, for the conhgurations of 

-h log M / d‘^X9/^A 


interest, we can take 39,40 


Set = j (^K - 6 -F ^e^'^dixdA - "0^ (^1 + 


2 logu 


(3.13) 


where 9* = and A is the conformal anomaly in the axion-dilaton-gravity system 
given by 40 


A^-e**\ax\' 

6 


(3.14) 
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Note that here we have expressed the X scalar held in terms of the canonically 

_t 

normalised scalar -0 dehned by X = e , which we will continue to use throughout 
this section. We note that the l/logw term is only relevant for solutions where the 


X-held is sourced, which are only briehy discussed in section 3.5 


The expectation value of the stress energy tensor is obtained by taking the func¬ 


tional derivative of the total action with respect to the boundary metric 41,42 


= lim 


( 2 K - 6 -t- ^e^^dixd^x - V’^(l + 


2 log M' 


+ logii (aj"’ - ^e‘*d‘xd>x{Sxf 


(3.15) 


where Kij is the extrinsic curvature of a m = constant hypersurface, and 5* = 

Using the boundary expansion of the helds in the previous section, we hnd the expec¬ 
tation value of the stress-energy tensor has the following non-vanishing components: 


rj^xx _ rpyy _ 
rj-^ZZ _ 


23 2945 4 4 2 

— oj A — —04 - d — —(^2) 

7 4032 7^ ^ 




-Tt 


13 


2731 


7 ^^^ 4032' 


f(^ 2 ) 


(3.16) 


This result is consistent with when X 2 —)■ 0. 

Similarly, we can calculate the one-point functions of the theory in order to hnd 
the vacuum expectation of the helds. For the scalar helds we hnd that expectation 
values and sources are given by 


(C’x)= 0 i X{o)=az, 

143 R 

(O^) = -—a^ + ^{I3, + 2Xi), 0(0) = 0 , 

{0^) = -V6X2, 0(0) =0. (3.17) 

We can now easily check that the Ward identities for the theory are satished. 
Firstly, diheomorphism invariance gives us the conservation of the stress-energy tensor 

-|- {Otp)Xj(j)(o) + (C^x)^iT(o) + (C^V’)Vj0(o) = 0, (3.18) 

which in our case is simply V*Tjj = 0, and is trivially satished. Similarly, the invari¬ 
ance of the theory under Weyl transformations leads to the conformal Ward anomaly 

7)* = -(4-A^)0(o)(O^)+^, (3.19) 
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where A is the conformal anomaly. From (3.16) we have T? = a^^/G and hence 


6 ’ 


(3.20) 


in agreement with a direct calculation of (3.14). 

By analytically continuing the time coordinate via t = —ir and demanding reg¬ 
ularity of the metric at u = Uh, we hnd that he Hawking temperature of the black 
holes is given by 

?V2| 

(3.21) 


B 

T = 


\:fu 


47r 


The entropy density of the black holes, s, can be obtained from the area of the black 
hole horizon and since we have set IGttG = 1 we have 

p~ 4 

s = 47r—(3.22) 

K 

We can calculate the free-energy density, w, by calculating the total on-shell Euclidean 
action via WV 0 I 3 = TItotai\os- In fact using the results of we can immediately 
obtain 


w = E — Ts , 


(3.23) 


where E = T** as well as the Smarr formula 


E-Ts = -T^ 


(3.24) 


As we explain in the appendix, these results can also be obtained by explicitly writing 
the bulk action as a total derivative in two different ways. 

Finally, we note that we can determine how various quantities transform under 


the scaling given in (3.9). For example, the free-energy transforms as 

w —)■ \~^w — A\~^ log A . 


(3.25) 


One can check that the Smarr formula is invariant under (3.9). 


3.4 Numerical construction of the black hole solutions 


We construct the black hole solutions by numerically solving the ODEs (3.3)-(3.6), 


subject to the boundary conditions given in (3.8),(3.10). Recall that for a hxed 


dynamical scale the black hole solutions can be parametrised by T/a. In practise we 
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T/a 

Figure 1: Plot showing the expectation value of the operator dual to the scalar 
held V’ (recall X = for the black hole solutions. The blue line corresponds to 

the solution of , while the red lines indicate the two new branches of solution, one 
that exists for T <Tc and the other which appears to exist only for T > Tc. 


set a = 1 and use a numerical shooting method in which we shoot from both near 
the black hole horizon and the holographic boundary and then match in the middle. 

We hnd that a new branch of black hole solutions appears at the critical temper¬ 
ature Tc/a ~ 1.8 X 10“^. In fact we hnd that there exist two distinct branches of 
black hole solutions carrying X hair, one that exists for T < Tc and the other that 
exists for T > Tc, as illustrated in hgures [T] and 

The solutions with T <Tc are the physically relevant solutions. In particular, we 
see from hgurej^that for T < Tc, where the black hole solutions of are unstable, 
this new branch of solutions has lower free energy. We thus conclude that there is 
a phase transition which, moreover, is a continuous phase transition. We emphasise 
that these hairy black holes are associated with a spontaneous phase transition since 
the boundary conditions we imposed for the held X corresponded to the dual operator 
(Pp acquiring an expectation value with no source. From the point of view of the 
D = 5 model this new phase does not appear to break any more symmetries than 
the background black holes. In particular, one can see from that potential in (2.6) 
does not have, for example, a Z 2 symmetry. However, after uplifting to type IIB, 
following an earlier discussion we know that when the '0-field acquires an expectation 
value then the 50(6) global symmetry is spontaneously broken to 50(4) x 50(2). 

The black hole solutions with T >Tc appear to be “exotic hairy black holes”. In 
particular, they only seem|^ to exist for T > Tc, in contrast to black holes associated 
with a first order transition which start existing for T > Tc and then turn around 
at some maximum temperature before continuing down to lower temperatures. We 


^We have checked that this is true up to T/a ~ 3. 
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(T/a)/ 10 - 


Figure 2: Plot showing the free energy of the black hole solutions, relative to the free 
energy at the critical temperature, Wc- The red line is the new branch of solution, 
while the blue line is the solution of [^. The left panel shows that the of branch 
of black holes that exist for T < has lower free energy, and hence is thermody¬ 
namically preferred. The right panel shows the black holes with T > Tc are not 
preferred. 


also observe from hgure that these black holes have higher free energy than the 
black holes of and hence are not thermodynamically preferred. We note that 
such exotic hairy black holes have appeared in other holographic constructions, both 


bottom up 28 and top-down 29 . 


3.5 Critical Exponents 

Having shown that there is a continuous phase transition at Tc, we now investigate the 
critical exponents of the transition. Somewhat surprisingly, we hnd that the phase 
transition does not have the same critical exponents as the majority of holographic 
phase transitions. 

The simplest critical exponent to calculate is /3, which is dehned by 

~ (Tc - T)^ (3.26) 

For our phase transition, from our numerics we hnd that /3 = 1, differing from the 
standard value = 1/2. There are several other important critical exponents for a 
phase transitiorj^ For example, the behaviour of the specihc heat for T < dehnes 
the exponent a via 


(Tc-T)-". 


(3.27) 


'‘For a discussion in the context of holography see 


44 
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This can be read off from the behaviour of the difference between the free energies 
of the two phases via Aw ~ (Tc — We hnd that in our transition a = — 1 in 

contrast to the standard value of a = 0. The remaining critical exponents are hxed 
by a, f3 using scaling relations. For example we have 

7 = 2-a-2/3-\ 5 = (2-a)/3-^-l, (3.28) 


where 7 , 5 are dehned by 


d{0^) 

dip 


(T, - Tf, i, ~ {0,,y 


(3.29) 


For our black holes we obtain 7 = 1, 5 = 2 in contrast to the standard results of 

7 = 1,5 = 3 . 

As a check that the scaling relations are indeed satisfied, we carried out a direct 
calculation of the exponent 6. To do this we constructed a more general class of black 
hole solutions with a source for the operator dual to ip. This required changing the 
boundary conditions in (3.8|) to allow for terms of the form logM. Having done this 


(which requires some effort), it is possible to see how the behaviour of (Op) needs to 
be changed as one switches on the source, while keeping the temperature fixed to be 
at the value T = Tc- Carrying out this procedure we found 5 = 2 in agreement with 
above. To summarise, the critical exponents for the new phase of black holes are 


(a,/ 3 , 7 , 5 ) = (-1,1,1,2), 


(3.30) 


in contrast to the standard values (a, 7 , 5) = (0,1/2,1, 3). Recall that the standard 

values arise from a Landau-Ginzburg model with quadratic and quartic terms in the 


free energy. In appendix B we discuss how the exponents (3.30) are associated with 


a free energy for a scalar order parameter with quadratic and cubic terms. 

Some bottom-up holographic superconducting phase transitions have been stud¬ 


ied in the probe approximation 30 32 and with back-reaction 33 , which also exhibit 


non-standard critical exponents. The specific critical exponents that we have found 
for our new black holes have also been found in models with a potential which con¬ 
tained terms that are cubic in the modulus of a complex scalar field. Although such 
couplings are rather unnatural for a charged scalar held, and it is difficult to see how 
they would arise form a top-down setting, we hnd that cubic terms in the potential 
for the neutral scalar held ip in our model are responsible for the non-mean held 
behaviour. 
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To see this we first note that if we expand our Lagrangian for the scalar held V’ 
around "0 = 0 we have 

+ 12 + 2ij‘^ - (3.31) 

In particular, the absence of a Z 2 symmetry -0 —'ip allows for the cubic term. We 
can contrast this model with a bottom up model in which the cubic term is absent: 

+ 12 + 2V’' + + O(V'), (3.32) 

with the rest of the Lagrangian unchanged. We have constructed the back-reacted 
black holes for this model and we hnd that the critical exponents for the phase 
transition now take the standard values. Furthermore, we have checked that varying 
the quartic terms does not change this result. 

3.6 Zero temperature scaling solution 

To investigate the low temperature behaviour of black hole solutions it is often illu¬ 
minating to examine the low temperature behaviour of Ts'/s, since if it approaches a 
constant it indicates an emergent scaling behaviour which one can then try to iden¬ 
tify. For the black hole solutions constructed in with X = 0 one hnds that s 
scales as s ~ tI and this is exactly the scaling behaviour that is associated with the 
Lifshitz-like anisotropic geometry found by that appears at T = 0 in the far IR. In 
fact this scaling behaviour is approximately present at the critical temperature phase 
transition as we see from hgure 

For our new black hole solutions with X 7 ^ 0 we also see from hgure that at very 
low temperatures s ~ T^. We therefore look for the existence of a scaling solution 
to the equations of motion of the form 

i3{u) = Bou^% X{u) = Xqm^^ . (3.33) 

However, by analysing the resulting algebraic equations one concludes that when 
X V 0 such solutions do not exist. Instead, we have found that the equations of 
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motion admit the following expansion as m —>■ oo: 

<A(n) _ _Jo _ 15232 




B{u) = 


X{u) = 


aM)4/9 10935 (aw) 16/90^9/2 

81 
m 
Bo 


4500 
3136i3o 


15/2 


(aM)2/3 3645(aM)20o^/^ 

4 28672 


+ .. 


+ ... , 
+ ... , 


+ 


(3.34) 


3(aM)4/90j/2 32805(aM)i6/9<^i4 

where Bo,4>o are constant. Furthermore, we have checked that the new black hole 
solutions start to approach this behaviour at low temperatures. Moreover, we can 
also show that this behaviour is associated with the observed scaling, s ~ Tt. 

To see this we hrst observe that the above expansion can be generalised to hnite 
temperatures, with the leading order expansion of X replaced with 

^ \ 28/9^ 

U \ 


Hu) = (^ ~ ( 


UhJ 


+ ..., 


(3.35) 


where Uh is the horizon radius. By combining (3.22) and the above hnite temperature 
solution, the entropy is given by 


s = 47ra^0Q ^^'^{auh) + 


while the Hawking temperature is given by 

r=^^^^(au,)“'/3 + .... 


IGtt 

Combining these two expressions we hnd, as claimed: 

4\^/^ 655367ri^/V 


^ 3 


218700 


49/4kj11/6 




+ .. 


(3.36) 


(3.37) 


(3.38) 


It is interesting to point out that the leading behaviour of the T = 0 solution 


given in (3.34) can be recast in the following form, after making the coordinate 
transformation u = ^ for some constant c: 

ds^ 


e7 ~ 


-2(3-9) , , 
p 3 \^dp 

X ~ p 


dp + dx'^ + dy'^ + p , 

-2/3 


(3.39) 


with 9 = —1, z = 2/?) and the bars denote that we have rescaled the coordinates. 
This is similar to the hyper-scaling solutions with Lifshitz exponent z and a hyper¬ 


scaling violation exponent 9 34 36 , but here the Lifshitz exponent is associated 
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Figure 3: Plots showing the temperature dependence of entropy, both along the high 
temperature branch of solution (right), and for the new low temperature branch of 
solution (left). The dotted lines are 11/3 and 8/3 on the left plot, and 3 and 8/3 for 
the right plot. We see that at very high temperatures s exhibits the standard scaling 
of the = 4 SYM plasma, before cooling towards a Lifshitz geometry. After the 
phase transition, the system moves away from the Lifshitz geometry, and towards the 
new scaling solution. 


with a spatial direction and not a time direction. Under the scaling {t, x, y, z, p) —)■ 
(At, Ah, A|/, A^/^A, Ap) we hnd that metric transforms as ds —>■ \ids. It is curious that 
the exponent z = 2/3 is the same value as for the unstable Lifshitz-like ground state 


with X = 1 given in (2.9). 


3.7 Thermal conductivity 

Having established the low temperature behaviour of the black holes, it is of interest 
to derive the DC thermal conductivity in the z-directioi]|^ k. To do so we follow the 


approach of 16 which showed how k can be obtained in terms of black hole horizon 
data. 


To make contact with it is convenient to write the black hole solutions in a 
slightly different form 


ds^ = —Udt^ + {dx^ + dy"^) + e^^dz '^, 

X = az, (j) = (j){r), X = X{r) 


(3.40) 


^Since the solutions are still translationally invariant in the x and y directions, the thermal 
conductivity in these directions is infinite. 
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where U, Vi and V 3 are functions of r. We assume that as r —>■ 00 , the functions have 
the following asymptotic form 

U^r^ + ..., + (3.41) 

and 0 —)■ 0 + o(m^), X —)■ 1 + o(m^). We now consider a small linearised perturbation 
about this class of black hole solutions that includes a piece that is linear in time: 


9tz = t 6 f 2 ir) + 5gta,^{r ), 

Qrz = e‘^^'^5hrz{r) , 

Xi = az + 6xi{r). 


(3.42) 


A key point is that this perturbation does not source the X-£eld. As a result the 


calculation of k is virtually unchanged from the derivation given in 16 . Rather than 
repeat the steps, we just quote the hnal result: 

dvrsT 


K = 


j2g2</> 


(3.43) 


J r=rh 


We showed in the previous section that the black holes with X 7 ^ 0 have s jnll/S 


at low temperatures. From (3.34) we can also determine the low temperature scaling 


behaviour of the dilaton to be [e^^)r=rf, ~ Hence the low temperature scaling 

of the thermal conductivity is given by k 2^10/3_ ggg |-pg PqIq 

solution is dual to a ground state that is thermally insulating in the direction of 
the linear axion. It is also interesting to contrast this result with the result for the 
(unstable) Lifshitz ground state, where k 


4 Discussion 

We have shown that the anisotropically deformed X = 4 Yang-Mills plasma studied 
in [^1^ has low temperature instabilities. The plasma undergoes a third-order phase 
transition, spontaneously breaking the global ^(^(G) symmetry down to S'0(4) x 
SO{2). We showed that critical exponents of the phase transition are given by 
(a, /3, 7 ,5) = (—1,1,1, 2) in contrast to the standard mean held theory values usually 
seen in holography. These critical values can be associated with a cubic Landau- 
Ginzburg free energy for a scalar order parameter, as discussed in appendix B. How¬ 
ever, such a free energy is unstable. In addition stabilising the free energy with a 
higher powers of the order parameter leads to a hrst order phase transition. By con¬ 
trast, in our holographic model the transition appears to be continuous, in fact third 
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order, provided that the branch of black holes with T > T^. does not turn around at 
some temperature and then go down to lower energies. Thus, our model underscores 
the difficulty in making a precise identihcation of the properties of the phase transi¬ 
tion just using a Landau-Ginzburg mean held approach. Perhaps a Landau-Ginzburg 
model with more helds might give a better description. It would be interesting to 
further clarify this point. 

It would also be interesting to know which critical exponents can be realised in 
string/M-theory constructions. In addition to the exponents that we found here it 
was shown that for a class of top-down i?-charged black holes the critical exponents 


are given by {a,/3,'y,6) = (1/2,1/2,1/2, 2) [^-47 . It is not clear if the bottom up 
constructions discussed in [30|-[33], which had more general exponents, can be embed¬ 


ded into top-down setting. Following 44 , it would also be of interest to explicitly 
calculate the dynamic critical exponents for our transition as well as others. 

We analysed the T = 0 limiting behaviour of the black holes describing the new low 
temperature phase. We showed that in the far IR there is an emergent leading order 
behaviour that is similar to the hyperscaling violation geometries but with spatial 
anisotropic scaling. This scaling behaviour implies that the thermal conductivity 
scales with temperature as k ~ yio/3 temperatures, revealing that the ground 

state is a thermal insulator. 

The black hole solutions were constructed using a consistent KK truncation that 
keeps a single scalar held X with = —2. This scalar held is part of a multiplet 
of twenty scalars that transform in the 20 ' of 50(6). All of these scalars become 
unstable at the critical temperature and it would be very interesting to investigate 
the full class of black hole solutions that emerges at Tc, which will generically break all 
of the 50(6) symmetry, and then follow them down to low temperatures. Although 
challenging, this could be investigated using the consistent truncation of that 
keeps twenty scalars parametrised by a symmetric, unimodular six by six matrix 


ij- 


As a hrst step one could analyse the truncation that keeps hve scalar helds, 
parametrised by the diagonal subset 49 , or even simpler, the truncations that keeps 


just two scalar helds 50 . The spontaneous breakdown of the global 50(6) symmetry 
will lead to Goldstone modes and it would also be interesting to study these further, 
both from the gravitational and held theory points of view. 

Finally, it is not difficult to show that the 0 = 5 model with metric, axion and 
dilaton (i.e. when X = 1) arises as a consistent truncation on an arbitrary hve- 
dimensional Sasaki-Einstein (SE) manifold, not just the hve-sphere. Therefore, the 
original black hole solutions of m also describe the high temperature phase of the 
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whole class of dual = 1 SOFT plasmas with an anisotropic deformation. For a 
given SE space, if there are no BF saturating modes in the spectrum then the black 
holes will not suffer the instabilities that we have described in this paper, and the 
Lifshitz ground state constructed in may be the true ground state of the system. 
On the other hand if there are BF saturating modes then the black holes will become 
unstable at some critical temperature. For a general SE manifold it is unlikely that 
there is a consistent truncation maintaining just one extra scalar held as we have 
studied in this paper for the case of the hve-sphere. This would mean that the 
corresponding black hole solutions describing the low temperature phase would need 
to be constructed directly in ten spacetime dimensions. Although this is likely to be a 
challenging task, it may be tractable to study the solutions near the phase transition 
and it would be particularly interesting to determine the critical exponents. 
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A Smarr relation 


We explain how to obtain the Smarr relation (3.24) via a direct calculation of the 


on-shell action. The bulk Euclidean bulk action is given by 

ruh 

hulk = -Arvols / duCbuik , 
where the Lagrangian density integrand is given by 
hbulk 


(A.l) 


= ^(^R- 3X-^{dXf + 4(X2 + 2X-1) - • (A.2) 

We would like to rewrite this as a total derivative in u after using the equations 
of motion. To achieve this we found it helpful to use the fact that after contraction 


of equation (2.2), we can write the integrand of the action as 






(A.3) 
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After some work we find that after using the equations of motion (3.3)-(3.6) the 
integrand can be written as 

^bulk 1 2 ^ h — — I , (A.4j 


U^\/B '2 m 3 


where the prime indicates differentiation with respect to the u coordinate. Using this 
expression we will get contributions to the on-shell action both from the horizon and 
the boundary. Using the near horizon and boundary expansions of the helds given 


in (3.8),(3.10), and combining with the boundary counter terms we deduce that the 


free energy density can be expressed as 


w = E — sT , 


(A.5) 


as in (3.23). 


On the other hand, using (3.3)-(3.6) we can also write the integrand in the form 

(A.6) 


Cbuik - I 2 - + — j 


This only gives contributions from the boundary leading to 

w = . 


(A.7) 


Combining these gives these expressions gives the Smarr relation (3.24). 


B 


Critical exponents for a cubic free ener 



Suppose we have a Landau-Ginzburg free energy functional for a scalar order param¬ 
eter, m, of the form 


^ ^ arn? brn? 

/ = /o+^ + ^ 


(B.l) 


with /o a constant, a = E, with t = {T — Tc)/Tc, and 6 is a temperature dependent 
constant which we take to be positive. We choose n so that a < 0 for T < E and 
we will be especially interested in the case n = 1. For the moment let us ignore the 
global instability for m < 0 and focus on the extrema at m = 0 and m = —a/b which 
exists when a < 0 i.e for T < T^. For the latter minimum we have m oc E and 
hence we conclude that /3 = n. To obtain a we want to differentiate the minimum 

®We would like to thank Makoto Natsuume for helpful discussions on this section. 
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value of the free energy with respect to T. Below we have / = /o + a^/ 66 ^ and 
hence we deduce that Td'^f/dT'^ oc ^ and thus a = 2 — 3n. Note that above 
Tc the free energy is constant and hence the specific heat vanishes. To determine 
S we add —mh to the free energy where h is a background source. We now have 
df /dm = am + — h and at T = T^., where a = 0, we deduce that the equilibrium 

configuration has m oc and hence 6 = 2. Finally, we consider the susceptibility 
X = dm/dh. At equilibrium we have am + bm?' — h = 0 and differentiating we deduce 
that X = l/(a + 26m). For T > T^. we have m = 0 and x = W"’, while for T < Tc we 
have m = —a/b and hence x = We thus deduce that 7 = n. When n = 1 the 

critical exponents are thus given by {a,(3,'j,6) = (— 1 , 1 , 1 , 2 ), exactly as we saw in 
our holographic phase transition. 

We now return to the issue of the global instability for m < 0. We first note that 
the instability would be eliminated if we were restricted to conhgurations with m > 0. 
Interestingly, the critical exponents that we have obtained were discussed in the 
context of a continuum generalisation of the Ashkin-Teller-Potts models associated 
with percolation problems, by imposing such a restriction [^. Note that we have 
no restrictions on the sign of the expectation value (Op), so this perspective is not 
available for our holographic phase transition. 

It is also worth pointing out that if we try to stabilise the free energy with higher 
powers of m, a quartic for example, then the model has a first order transition, again 
unlike what we see in our holographic transition. More explicitly we can add a term 
cm^/4 to the free energy in (B.l) with c > 0. Now for high temperatures, a > 6 ^/(4c), 
the free energy has a minimum at m = 0. For 26^/(9c) < a < b‘^/{4:c) there is an 
additional minimum at m = mi = —b/{2c) — [ 6 ^ — 4ac]^/^/( 2 c), which, has higher free 
energy than the minimum at m = 0. For 0 < a < 26^/(9c) the minimum at mi has 
lower free energy than the minimum at m = 0 and there is a first order transition at 
a = 26^/(9c). For a < 0, m = 0 becomes a maximum of the free energy with a new 
minimum appearing at m = m 2 = —b/{2c) + [ 6 ^ — 4ac]^/^/(2c). This m 2 minimum 
is the one associated with the critical exponents for the cubic with c = 0 that we 
discussed above, but it is simple to see that the mi minimum is always preferred. 

In summary, we see that while the cubic Landau-Ginzburg model for a single 
scalar order parameter in a certain sense gives rise to the critical exponents we see in 
our holographic phase transition it does not capture key features. Perhaps a model 
containing more fields might be more effective. 
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